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Flameless  combustion  is  able  to  provide  high  combustion  efficiency  with  low  NOx  and  soot  emissions. 
The  present  work  aims  at  investigating  the  role  of  closure  sub-models  for  the  modelling  of  a  flameless 
furnace,  as  well  as  the  main  NO  formation  paths.  Among  the  different  turbulence  models  that  were 
tested,  modified  k—e  provides  the  best  agreement  with  the  experimental  data,  especially  for  temperature 
measurements.  Reynolds  stress  model  leads  to  smaller  deviation  for  radial  velocity  predictions.  Since  in 
flameless  combustion  regime  the  turbulence-chemistry  interaction  as  well  as  the  kinetic  mechanism 
play  a  fundamental  role,  the  Eddy  Dissipation  Concept  (EDC),  coupled  with  four  different  kinetic  schemes 
(JL,  KEE58,  GRI  2.11  and  GR1  3.0)  was  considered.  The  GR1  2.11  and  KEE58  mechanisms  perform  better, 
thus  confirming  the  necessity  of  turbulence/chemistry  interaction  models  accounting  for  finite-rate 
chemistry  when  flameless  combustion  is  studied.  As  far  as  NO  emissions  are  concerned,  the  N20  in¬ 
termediate  NO  mechanism  is  found  to  play  a  major  role,  while  thermal  NO  formation  mechanism  is  not 
as  relevant  as  in  traditional  combustion  regime. 

An  assessment  of  the  uncertainty  related  to  the  choice  of  boundary  conditions  as  well  as  to  the  choice 
of  the  parameters  of  the  physical  models  is  also  performed.  Finally  the  operation  characteristics  (such  as 
the  recirculation  rate  and  the  location  of  the  reaction  zone)  of  the  furnace  are  evaluated. 

©  2014  Elsevier  Ltd.  All  rights  reserved. 


1.  Introduction 

Flameless  combustion  [1],  also  known  as  Moderate  and  Intense 
Low-Oxygen  Dilution  (MILD)  [2]  or  HiTAC  combustion  [3]  is  able  to 


*  Corresponding  author.  Universite  Libre  de  Bruxelles,  Service  d’Aero-Thermo- 
Mecanique,  Avenue  F.  D.  Roosevelt  50,  B-1050  Bruxelles,  Belgium.  Tel.:  +32  2  650 
26  80;  fax:  +32  2  650  27  10. 

**  Corresponding  author.  Universita  di  Pisa,  Dipartimento  di  Ingegneria  Civile  e 
Industriale,  Largo  L.  Lazzarino  2,  1-56122  Pisa,  Italy.  Tel.:  +39  050  22  17  897; 
fax:  +39  050  22  17  866. 

E-mail  addresses:  chiara.galletti@unipi.it  (C.  Galletti),  Alessandro.Parente@ulb. 
ac.be  (A.  Parente). 

http://dx.doi.Org/10.1016/j.applthermaleng.2014.ll.005 

1359-4311/©  2014  Elsevier  Ltd.  All  rights  reserved. 


provide  high  combustion  efficiency  with  low  NO*  and  soot  emis¬ 
sions.  The  increasing  interest  in  flameless  combustion  is  motivated 
by  the  large  fuel  flexibility,  representing  a  promising  technology  for 
low-calorific  value  fuels  [4],  high-calorific  industrial  wastes  as  well 
as  in  presence  of  hydrogen  5  .  Moreover,  flameless  combustion  is 
very  stable  and  noiseless  [6],  so  it  could  be  potentially  suited  for  gas 
turbine  applications  [7]  where  conventional  operations  may  lead  to 
significant  thermo-acoustic  instabilities  (“humming”)  and  stresses. 

Flameless  combustion  needs  the  reactants  to  be  preheated 
above  their  self-ignition  temperature  and  enough  inert  combustion 
products  to  be  entrained  in  the  reaction  region,  in  order  to  dilute 
the  flame.  As  a  result,  the  temperature  field  is  more  uniform  than  in 
traditional  non-premixed  combustion  systems,  and  it  does  not 
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show  high  temperature  peaks.  Hence,  NO*  formation  is  suppressed 
as  well  as  soot  formation,  due  to  the  lean  conditions,  low  temper¬ 
atures  and  the  large  CO2  concentration  in  the  exhausts. 

This  combustion  regime  appears  to  be  still  worthy  of  further 
investigations  and  attention.  In  particular,  the  fundamental  mech¬ 
anism  of  the  interaction  between  turbulent  mixing  and  chemical 
kinetics  needs  to  be  elucidated.  With  respect  to  conventional 
flames,  turbulence  levels  are  enhanced  (due  to  the  high  momentum 
of  the  gases),  thus  mixing  timescales  are  reduced;  on  the  contrary 
chemical  timescales  are  increased  due  to  dilution  of  the  reactants 
[8],  In  flameless  combustion,  the  Damkohler  number  approaches 
unity  [9],  implying  that  both  mixing  and  chemical  kinetics  should 
be  taken  into  account  when  modelling  such  a  regime,  resulting  in  a 
very  challenging  problem.  In  addition,  most  of  the  available  models 
have  been  derived  for  conventional  flames;  hence  they  need  to  be 
validated  and  eventually  revised  for  non-conventional  regimes. 

From  a  computational  perspective,  the  role  of  the  combustion 
model  and  the  possible  degree  of  simplification  of  chemical  kinetics 
have  not  been  rigorously  and  systematically  assessed  in  the  context 
of  flameless  combustion.  Encouraging  results  in  literature  on  the 
modelling  of  such  flames  have  been  found  using  the  Eddy  Dissi¬ 
pation  Concept  (EDC)  model  by  Magnussen  [10],  coupled  with 
detailed  chemical  mechanism.  Such  modelling  strategy  has  been 
successfully  applied  to  the  Jet  in  Hot  Coflow  burner  by  Christo  and 
Dally  [11  ]  and  Aminian  et  al.  [12],  The  same  approach  leads  to  good 
results  also  in  presence  of  hydrogen  in  the  fuel,  as  reported  by 
Parente  et  al.  [13]  and  Galletti  et  al.  [14  for  a  self-recuperative 
burner  and  a  lab-scale  burner,  respectively,  operating  in  flameless 
combustion  conditions. 

However,  recent  investigations  carried  out  by  De  et  al.  [15]  and 
Aminian  et  al.  [16]  have  indicated  the  need  for  revising  the  EDC 
constants  in  the  framework  of  flameless  combustion  modelling. 

All  the  works  reported  above  show  that  the  oxidation  scheme 
may  strongly  impact  the  results,  as  indicated  by  Shabanian  et  al. 

[17] .  The  global  schemes  are  unsuited  and  generally  lead  to  a  strong 
over  prediction  of  the  flame  temperatures.  Recently  Parente  et  al. 

[18]  applied  Principal  Component  Analysis  to  a  set  of  measure¬ 
ments  taken  in  flameless  combustion  conditions,  showing  that 
approaches  based  on  single  progress  variable  are  not  suited  for  the 
description  of  such  combustion  regime  and  finite-rate  chemistry 
models  are  needed. 

Finally,  the  modelling  of  NO  emissions  deserves  also  special 
attention.  The  Zeldovich  thermal  mechanism  is  the  major 
contributor  to  NO  in  most  of  the  conventional  combustion  system; 
however,  in  flameless  combustion,  the  lower  temperatures  and  the 
absence  of  large  fluctuations  inhibit  NO  formation  through  such  a 
mechanism.  As  a  result,  NO  emissions  are  controlled  by  other  for¬ 
mation  routes,  such  as  the  Fenimore's  prompt  NO  and/or  N2O  in¬ 
termediate  [19],  Therefore,  it  is  necessary  to  incorporate  all 
potentially  relevant  formation  paths  in  the  numerical  model. 
Moreover,  other  routes  may  become  relevant  with  non- 
conventional  fuels,  such  as  the  NNH  pathway  for  H2  containing 
fuels  [14],  The  prediction  of  NO  formation  in  flameless  combustion, 
at  low  temperatures  and  with  high  concentration  of  H2  in  the  fuel 
stream  has  been  studied  by  Parente  et  al.  [20],  They  found  that  the 
inclusion  of  non-conventional  NO  formation  routes,  i.e.  N20  inter¬ 
mediate  and  NNH,  is  crucial  for  characterizing  the  pollutant 
emissions. 

The  elucidation  of  the  above  topics  needs  high  fidelity  and 
comprehensive  experimental  data  to  validate  the  numerical 
models.  The  Jet  in  Hot  Coflow  (JHC)  burner  [21  ,  the  Delft  Jet  in  Hot 
Coflow  (DJHC)  [22,23]  and  the  Cabra  flame  [24]  have  been 
conceived  to  emulate  flameless  conditions  by  feeding  diluted  and 
hot  streams  to  the  burner.  They  constitute  a  strong  asset  for  the 
validation  of  numerical  models  as  they  have  been  equipped  with 


advanced  diagnostics  to  measure  mean  and  fluctuating  variables 
(e.g.  chemical  species,  temperature,  velocities). 

However,  in  the  industrial  practice,  flameless  conditions  are 
obtained  by  means  of  the  massive  internal  recirculation  of  flue 
gases,  which  allows  diluting  the  fresh  gases  before  they  reach  the 
reaction  zone.  Such  recirculation  is  generally  achieved  through 
special  designs  of  the  feeding  jets  as  well  as  of  the  combustion 
chamber.  A  recent  review  of  different  designs  of  flameless  com¬ 
bustors  is  provided  by  Arghode  and  Gupta  [25,26], 

The  recirculation  affects  both  mixing  and  chemical  timescales, 
so  that  conceptually  these  burners  are  different  from  JHC,  DJHC  and 
Cabra  flames,  which  act  solely  on  the  chemical  timescale.  A  few 
experimental  investigations  of  flameless  furnaces,  based  on  inter¬ 
nal  recirculation  of  exhaust  gases  can  be  found  in  literature.  Szego 
et  al.  [19]  described  the  performance  and  stability  characteristics  of 
a  parallel  jet  flameless  combustion  burner  system  in  a  20  kW 
laboratory-scale  furnace.  Mi  et  al.  [27]  investigated  a  20  kW  recu¬ 
perative  MILD  furnace,  using  EDC  combustion  model  with  global 
kinetic  schemes  for  methane  and  ethane.  Plessing  et  al.  [28]  and 
Ozdemir  and  Peters  [29]  provided  a  useful  set  of  experimental  data 
(velocity,  temperature  and  NO  emissions)  on  a  5.4  kW  furnace  fed 
with  methane  and  operating  in  flameless  regime,  subsequently 
modelled  by  Coelho  and  Peters  [30]  using  the  Eulerian  Particle 
Flamelet  model.  Their  results  showed  some  discrepancies  in  the 
prediction  of  flow  field  as  well  as  the  overestimation  of  NO  levels  at 
the  outlet  section.  Dally  et  al.  [31  extended  the  investigation  of  the 
same  furnace  to  more  fuels  and  equivalence  ratios.  Verissimo  et  al. 
[32]  experimentally  investigated  a  10  kW  reversed-flow  cylindrical 
furnace,  for  which  simulation  were  performed  by  Graqa  et  al.  [33], 
The  Authors  compared  the  EDC  model  coupled  with  the  DRM-19 
mechanism  and  the  composition  PDF  (C-PDF)  model,  showing  a 
general  good  agreement  between  predictions  and  experiments. 
Danon  et  al.  [34]  performed  a  parametric  study  on  a  300  l<Wth 
furnace  equipped  with  three  pairs  of  regenerative  flameless  com¬ 
bustion  burners  with  the  objective  of  optimizing  the  furnace  per¬ 
formance.  These  experimental  results  were  used  as  validation  data 
for  a  set  of  Computational  Fluid  Dynamics  (CFD)  simulations  of  the 
furnace  reported  in  Ref.  [35j.  The  authors  showed  that  the  EDC 
model  in  combination  with  the  realizable  k—e  model  and  a  skeletal 
chemistry  mechanism  allowed  reproducing  the  main  furnace  per¬ 
formance  for  all  the  investigated  burner  configurations.  Rebola 
et  al.  [36]  performed  an  experimental  investigation  on  a  small-scale 
flameless  combustor,  defining  the  range  of  operating  conditions 
allowing  operating  in  flameless  conditions.  Cameretti  et  al.  [37] 
discussed  some  aspects  related  to  the  employment  of  liquid  and 
gaseous  bio-fuels  in  a  micro-gas  turbine  operating  in  flameless 
regime,  showing  numerically  the  energetic  and  environmental 
advantages  related  to  the  use  of  those  fuels.  Recently  Huang  et  al. 
[38]  studied  the  emissions  from  a  flameless  combustion  staged 
combustor.  The  authors  found  that  the  flameless  regime  yields 
lower  NO  emissions  compared  to  the  traditional  diffusion  com¬ 
bustion  mode,  and  the  N2O  intermediate  mechanism  dominates 
the  NO  production. 

The  present  work  aims  at  investigating  the  role  of  closure  sub¬ 
models  for  the  modelling  of  the  flameless  furnace  of  Plessing 
et  al.  [28]  and  Ozdemir  and  Peters  [29],  as  well  as  at  identifying  the 
main  NO  formation  paths.  The  selected  experiment  set-up  was 
chosen  as  it  represents  an  optimal  compromise  between  lab-scale 
and  industrial  systems.  It  shows,  in  fact,  the  typical  feature  of  in¬ 
dustrial  flames  systems,  i.e.  the  internal  aerodynamic  recirculation, 
allowing,  at  the  same  time,  a  sufficiently  detailed  characterization 
of  the  system  performances. 

While  existing  literature  has  pointed  out  the  crucial  role  of  finite 
rate  chemistry  models  and  detailed  kinetics  in  flameless  regime, 
little  emphasis  has  been  devoted  to  the  quantification  of  the 
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uncertainty  related  to  the  boundary  conditions  and  physical 
models.  The  analysis  focuses  on  the  quantitative  assessment  of  the 
scenario,  i.e.  boundary  conditions,  and  modelling  uncertainties  on 
the  results.  In  conjunction  with  appropriate  validation  metrics,  this 
allows  identifying  the  most  sensitive  parameters  for  the  numerical 
simulations  and  developing  a  predictive  model,  which  is  accurate 
enough  for  the  description  of  the  system.  Finally  the  operation 
characteristics  (such  as  the  recirculation  rate  and  the  location  of  the 
reaction  zone)  of  the  actual  furnace  are  evaluated. 

2.  Test  case 

A  sketch  of  the  combustion  chamber  is  provided  for  sake  of 
clarity  in  Fig.  1,  even  though  all  details  of  the  system  and  experi¬ 
mental  campaigns  can  be  found  in  Refs.  [28,29].  The  combustion 
chamber  is  a  parallelepiped  (with  250  x  250  mm2  cross-section 
and  485  mm  length)  with  both  the  burner  and  the  exhaust  outlet 
placed  at  the  bottom.  The  multi-nozzle  burner  consists  of  a  central 
nozzle  (inner  diameter,  i.d.  =  4.7  mm)  for  fuel  inlet;  this  nozzle  is 
conically  elevated  25  mm  from  the  6  peripheral  nozzles 
(i.d.  =  5  mm)  for  air  inlet,  located  40  mm  away  from  the  centre.  The 
air  nozzles  are  located  16  mm  higher  than  the  15.5  mm  wide 
annular  exit  (i.d.  =  93  mm)  for  exhaust  gases.  The  fuel  is  methane. 
Measurements  were  made  for  air  and  fuel  mass  flow  rates  of  6.5 
and  0.38  kg/h,  respectively,  which  correspond  to  an  overall  equiv¬ 
alence  ratio  cp  =  1.  Both  fuel  stream  and  air  stream  were  preheated, 
to  650  I<  and  1150  K  respectively. 

3.  Numerical  simulations 

The  CFD  model  was  defined  using  the  commercial  code  Fluent 
6.3  by  Ansys  Inc.  The  geometry  and  grid  were  realized  using  the 
software  Gambit.  Due  to  the  combustion  chamber  symmetry,  just  a 
quarter  of  the  geometry  is  considered  for  the  numerical  modelling. 
The  grid,  chosen  after  a  mesh  independency  study,  contained 
501,000  cells;  tetrahedrons  were  used  near  the  burner  zone, 
whereas  hexahedral  cells  were  used  in  the  remainder  of  the 
furnace.  The  choice  of  a  hybrid  grid  was  aimed  at  reducing  the 
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Fig.  1.  Sketch  of  the  experimental  furnace. 


number  of  cells  required  to  discretize  the  computational  domain 
accurately. 

3.1.  Boundary  conditions 

As  far  as  boundary  conditions  are  concerned,  mass  flow  condi¬ 
tions  were  specified  at  both  air  inlets  and  fuel  inlets;  air  inlet 
temperature  is  1150  K  whereas  fuel  inlet  temperature  is  650  K.  A 
pressure  outlet  condition  was  employed  for  the  exhaust  gases  exit. 
As  far  as  walls  are  concerned,  they  were  considered  isothermal.  The 
temperature  of  1313  K  was  chosen,  as  suggested  in  Ref.  29],  For 
flameless  condition  operations,  this  appears  an  acceptable 
approximation,  as  supported  in  Ref.  1  ].  Simulations  with  fixed  wall 
temperature  were  preferred  over  imposed  heat  flux  boundary 
condition,  since  the  latter  exhibited  numerical  convergence  prob¬ 
lems,  leading  to  unstable  solutions  and  to  the  extinction  of  the 
combustion  process.  Simulations  were  also  run  considering  a 
convective  heat  transfer  at  the  wall,  in  order  to  verify  the  goodness 
of  the  isothermal  hypothesis.  The  mean  heat  transfer  coefficient, 
calculated  from  a  global  heat  balance  on  the  furnace,  was  evaluated 
at  7  W/m2  K,  which  is  compatible  with  a  natural  convection 
process. 

3.2.  Physical  models 

Favre-averaged  Navier— Stokes  equations  were  solved,  using  the 
following  turbulence  models  to  resolve  Reynolds  stresses: 

-  Standard  k—e  turbulence  model  (SKE)  [39]; 

-  Modified  k—e  model  (MKE)  proposed  by  Morse  [40]  consisting 
in  the  variation  of  the  first  constant  of  the  turbulent  kinetic 
energy  dissipation  equation  from  1.44  to  1.6; 

-  Renormalization  group  (RNG)  k—e  [41]; 

-  Realizable  k-e  (RKE)  [42]; 

-  Reynolds  Stress  (RSM)  model  [43], 

Radiation  effects  were  accounted  for  using  the  Discrete  Ordinate 
(DO)  radiation  model  with  the  Weighted  Sum  of  Gray  Gases 
(WSGG)  model  for  the  participating  media  radiation,  using  the 
coefficients  proposed  by  Smith  et  al.  [44], 

Turbulence— chemistry  interactions  were  modelled  using  EDC 
[10].  As  mentioned  in  the  introduction,  the  EDC  model  is  able  to 
account  for  finite-rate  effects  and  thus  it  can  incorporate  detailed 
kinetic  schemes.  Within  the  present  study,  four  different  kinetic 
schemes  were  used  for  methane  oxidation: 

-  JL  scheme,  which  is  a  4-step  global  kinetic  mechanism  of  Jones 
and  Lindstedt  [45],  used  both  for  methane  and  hydrogen- 
enriched  fuels.  It  involves  7  species; 

-  KEE58  scheme,  which  is  a  skeletal  mechanism  made  of  17  species 
and  58  chemical  reversible  reactions  [46]; 

-  GR1  2.11  scheme,  which  is  a  detailed  mechanism.  It  was  imple¬ 
mented  without  the  NO*  reactions,  resulting  in  175  chemical 
reactions  involving  31  species  [47]; 

-  GRI  3.0  scheme,  which  is  a  detailed  mechanism.  It  was  imple¬ 
mented  without  the  NO*  reactions,  resulting  in  217  chemical 
reactions  involving  35  species  [48], 

The  In-Situ  Adaptive  Table  (ISAT)  [49]  was  coupled  to  EDC  to 
reduce  the  computational  costs.  An  error  tolerance  of  10  5  was 
selected  to  obtain  table-independent  results. 

As  for  NO  emissions,  different  formation  routes  were  consid¬ 
ered:  thermal  NO  mechanism,  Fenimore's  Prompt  mechanism  and 
intermediate  N2O  mechanisms.  The  thermal  NO  formation  is 
modelled  using  a  Finite  Rate  combustion  model  with  a  simplified 
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one-step  mechanism,  obtained  from  the  Zeldovich  scheme  by 
assuming  a  steady  state  for  the  N  radicals  and  relating  the  O  radical 
concentration  to  that  of  oxygen  by  means  of  the  dissociation  re¬ 
action  50],  The  prompt  NO  formation  is  modelled  following  De 
Soete  [51],  Finally,  the  N20  intermediate  mechanism  is  modelled 
according  to  the  equation  proposed  by  Malte  and  Pratt  [52], 

All  the  above  NO  formation  kinetic  rates  are  integrated  over  a 
probability  density  function  (PDF)  for  temperature,  to  account  for 
the  effect  of  temperature  fluctuations  on  the  mean  reaction  rates. 
The  assumed  PDF  shape  is  that  of  a  beta  function  [53]  and  is  eval¬ 
uated  through  the  temperature  variance,  the  latter  solved  by  means 
of  a  transport  equation. 


the  evaluation  of  the  intrinsic  experimental  error  and  the  aleatory 
nature  of  the  phenomenon  under  evaluation.  The  experimental 
data  provided  in  Ref.  [29]  do  not  allow  quantifying  these  quantities 
in  details,  especially  the  statistical  variability  of  the  measurements, 
which  requires  a  collection  of  multiple  observations.  Therefore, 
information  from  the  literature  [55]  were  used  to  estimate  poten¬ 
tial  uncertainties  in  experimental  observations,  i.e.  wall  tempera¬ 
ture,  to  evaluate  the  effect  of  uncertainties  in  boundary  conditions 
on  the  results,  as  reported  in  Section  5. 

4.  Results 

4.3.  Effect  of  turbulence  model 


3.3.  Numerical  settings 

A  second-order  upwind  discretization  scheme  was  used  for  all 
equations  and  the  SIMPLE  algorithm  was  employed  for  pressur¬ 
e-velocity  coupling.  The  simulations  were  run  until  the  residuals 
for  all  the  resolved  quantities  levelled  out,  resulting  in  a  decrease  of 
at  least  six  orders  of  magnitude.  In  addition  flow  field  variables  at 
different  locations  were  monitored  to  check  convergence  to  the 
steady  state  solution. 


3.4.  Validation  metrics 


In  order  to  quantitatively  measure  the  agreement  between 
experimental  data  and  computational  predictions,  validation  met¬ 
rics  were  employed,  following  Oberkampf  and  Barone  [54], 

For  temperature,  the  average  error  metric  normalized  with 
respect  to  the  measured  value  was  evaluated,  which  is  defined  as: 
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where  Ye  and  Ym  are  the  mean  measurement  and  the  predicted 
value  of  variable  Y. 

For  velocity  measurements,  it  was  chosen  to  refer  to  an  absolute 
average  error  metric,  defined  as: 
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The  choice  of  an  absolute  error  metric  for  velocity  measure¬ 
ments  was  justified  by  the  need  of  avoiding  division  by  zero  in  the 
error  calculation.  However,  this  kind  of  metric  may  not  be  enough 
to  assess  the  goodness  of  the  several  models  used.  To  better  un¬ 
derstand  the  deviation  of  predictions  from  measurements  the 
Normalized  Root-Mean-Square  Error  ( nrmse )  was  also  computed. 
This  metric  represents  the  sample  standard  deviation  of  the  dif¬ 
ferences  between  the  value  predicted  by  a  model  and  the  value 
actually  observed,  divided  by  the  range  of  observed  values  of  the 
variable  being  predicted  and  it  is  defined  as: 


V  N 

nrmse  =  = - = -  (3) 

Ye, Max  ~  Ye jyijn 

In  addition,  the  quantitative  agreement  between  predictions 
and  experiments  was  also  assessed  with  the  coefficient  of  deter¬ 
mination  (R2),  from  the  parity  plot  of  the  axial  velocity  and  the 
velocity  fluctuations. 

Beside  the  metrics  provided  by  Equations  (1)— (3),  the  con¬ 
struction  of  an  error  validation  metric  requires  the  estimation  of 
the  interval  containing  the  true  error.  The  latter  is  obtained  from 


Figs.  2  and  3  show  the  comparison  between  experimental  radial 
profiles  of  axial  and  radial  velocity,  respectively,  and  those  pre¬ 
dicted  at  different  axial  coordinates  with  different  turbulence 
models  and  using  the  EDC  model  with  KEE58  kinetic  scheme. 

The  radial  profile  of  axial  velocity  obtained  very  close  to  the 
burner  exit  (z  =  0.012  m)  shows  an  excellent  agreement  with 
experimental  data  with  all  turbulence  models.  At  further  distances 
the  predicted  data  capture  the  trend,  but  typically  predicted  pro¬ 
files  are  shifted  towards  larger  distances.  For  the  radial  velocity, 
predicted  values  are  generally  lower  that  the  measured  ones  for 
z  =  0.012,  0.112  and  0.212  m,  although  the  trend  of  the  profile  is 
captured. 

Tables  1  and  2  list  the  absolute  average  validation  metrics  as 
well  as  the  nrmse  for  the  profiles  of  axial  and  radial  velocities  as 
obtained  with  the  different  turbulence  models.  None  of  the  tur¬ 
bulence  models  tested  gives  perfect  agreement  with  experimental 
data.  Moreover  the  behaviour  of  the  different  model  is  similar,  with 
global  absolute  deviation  from  the  measured  value  ranging  from 
2.84  to  4.06  m/s  for  the  axial  velocity  (minimum  value  being 
0.58  m/s  with  MKE  and  maximum  value  being  10.17  m/s  with  RNG) 
and  from  0.92  to  0.99  m/s  for  the  radial  velocities  (minimum  value 
being  0.02  m/s  with  SKE  and  maximum  value  being  2.37  m/s  with 
RNG).  The  same  conclusions  may  be  drawn  considering  the  nrmse. 
The  models  that  perform  better  are  the  SKE  (0.31  for  axial  velocity 
and  0.59  for  the  radial  one)  and  the  RSM  (0.29  for  axial  velocity  and 
0.61  for  the  radial  one). 

The  comparison  between  experimental  radial  profiles  of  the  axial 
component  of  velocity  fluctuations  and  those  predicted  with  different 
turbulence  models  is  reported  in  Fig.  4.  Standard  k—e  and  RSM  show 
the  best  agreement  with  the  experimental  data,  even  if  some 
discrepancy  can  be  noted.  The  absolute  average  validation  metrics  as 
well  as  the  nrmse  are  shown  in  Table  3.  Standard  k—e  is  indeed  the 
model  that  provides  the  lowest  deviation  (0.5  m/s  and  0.41  for  nrmse), 
even  though  the  other  models  show  a  quite  similar  behaviour  (de¬ 
viations  ranging  from  0.5  to  0.87  m/s,  minimum  value  being  0.01  m/s 
with  SKE  and  maximum  value  being  4.41  m/s  with  MKE). 

To  confirm  these  observations,  the  coefficient  of  determination 
(R2)  for  the  parity  plot  of  the  measured  and  calculated  velocities 
(Fig.  5)  is  reported  in  Table  4.  The  standard  k—e  model  shows  the 
best  agreement  with  the  experimental  results  (R2  =  0.893  for  axial 
mean  velocity  and  R2  =  0.869  for  its  fluctuation).  As  far  as  the  radial 
velocity  is  concerned,  the  data  show  a  low  degree  of  correlation, 
also  imputed  to  the  small  absolute  values  of  velocity. 

The  prediction  of  the  velocity  field  presented  here  is  in  line  with 
the  results  provided  by  Coelho  and  Peters  [30  ,  who  reported  dis¬ 
crepancies  between  the  measured  and  predicted  velocity  profiles. 
However,  the  lack  of  a  discussion  of  the  intrinsic  uncertainties 
related  to  velocity  measurements  does  not  allow  completely 
assessing  the  performances  of  the  different  turbulence  models  by 
means  of  the  analysis  of  the  predicted  and  measured  velocity  fields. 
Additional  information  is  needed  to  select  the  turbulence  model  for 
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Fig.  2.  Comparison  between  experimental  axial  velocity  profiles  and  those  predicted  with  different  turbulence  models,  at  different  axial  coordinates.  EDC  combustion  model,  KEE58 
kinetic  mechanism. 


the  subsequent  investigation  of  the  effect  of  combustion  models 
and  kinetic  mechanisms. 

Fig.  6(a)  shows  the  measured  temperature  profiles  along  the 
combustion  chamber  axis  and  those  obtained  using  the  different 
turbulence  models  (with  EDC  and  KEE58).  The  modified  k— e  model 
provides  the  best  results.  Standard  k—e,  RKE  and  RSM  underesti¬ 
mate  the  temperature  along  the  axis,  while  RNG  yields  higher 
temperature.  The  discrepancies  in  the  temperature  profile  along 


Radial  Distance  [m] 


the  axis,  provided  by  the  different  models  are  due  to  a  completely 
different  fluid  dynamics  of  the  flame.  Relative  error  metrics  were 
used  in  this  case  to  quantitatively  evaluate  the  agreement  between 
experimental  and  predicted  data.  The  modified  k—e  model  leads  to 
an  average  deviation  of  about  2%,  thus  lower  than  the  other  models 
(6%  for  RNG,  7%  for  RKE  and  8%  for  SKE  and  RSM).  Based  on  these 
considerations,  the  MKE  model  was  chosen  for  the  following 
analysis. 


Radial  Distance  [m] 


Radial  Distance  [m] 


Fig.  3.  Comparison  between  experimental  radial  velocity  profiles  and  those  predicted  with  different  turbulence  models,  at  different  axial  coordinates.  EDC  combustion  model,  KEE- 
58  kinetic  mechanism. 
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Table  1 

Average  absolute  validation  metrics  [m/s]  and  normalized  root-mean-squared  error 
[— ]  for  axial  velocity  with  different  turbulence  models.  EDC,  KEE58. 

z  [m]  Average  absolute  metrics  [m/s]  Normalized  root-mean-squared  error  [-] 


MKE 

SKE 

RNG 

RKE 

RSM 

MKE 

SKE 

RNG 

RKE 

RSM 

0.012 

2.77 

2.25 

2.99 

2.82 

2.18 

0.09 

0.08 

0.10 

0.09 

0.08 

0.112 

4.56 

4.64 

5.90 

4.27 

4.20 

0.45 

0.40 

0.51 

0.39 

0.37 

0.212 

3.95 

3.23 

3.25 

4.67 

2.82 

0.40 

0.30 

0.40 

0.36 

0.27 

0.312 

5.00 

2.79 

3.19 

4.41 

2.60 

0.87 

0.47 

0.53 

0.77 

0.44 

Total 

4.06 

2.84 

3.42 

3.95 

2.95 

0.45 

0.31 

0.38 

0.40 

0.29 

Table  3 

Average  absolute  validation  metrics  [m/s]  and  normalized  root-mean-squared  error 
[— ]  for  axial  component  of  velocity  fluctuations  with  different  turbulence  models. 
EDC,  KEE58. 


z  [m]  Average  absolute  metrics  [m/s]  Normalized  root-mean-squared  error  [— ] 


MKE 

SKE 

RNG 

RKE 

RSM 

MKE 

SKE 

RNG 

RKE 

RSM 

0.012 

1.40 

0.90 

1.35 

1.32 

1.38 

0.25 

0.15 

0.22 

0.22 

0.23 

0.112 

0.46 

0.55 

0.73 

0.63 

0.55 

0.25 

0.30 

0.42 

0.34 

0.29 

0.212 

0.90 

0.54 

0.33 

1.03 

0.62 

0.67 

0.45 

0.25 

0.80 

0.51 

0.312 

0.77 

0.40 

0.15 

0.82 

0.42 

1.32 

0.72 

0.36 

1.49 

0.73 

Total 

0.77 

0.50 

0.62 

0.80 

0.61 

0.62 

0.41 

0.31 

0.71 

0.44 

Table  2 

Average  absolute  validation  metrics  [m/s]  and  normalized  root-mean-squared  error 
[— ]  for  radial  velocity  with  different  turbulence  models.  EDC,  KEE58. 


z  [m]  Average  absolute  metrics  [m/s]  Normalized  root-mean-squared  error  [-] 


MKE 

SKE 

RNG 

RKE 

RSM 

MKE 

SKE 

RNG 

RKE 

RSM 

0.012  1.31 

1.48 

1.16 

1.38 

1.08 

0.67 

0.78 

0.59 

0.70 

0.55 

0.112  0.39 

0.23 

0.31 

0.31 

0.23 

0.71 

0.42 

0.57 

0.56 

0.47 

0.212  0.53 

0.50 

0.51 

0.49 

0.60 

0.69 

0.68 

0.68 

0.64 

0.81 

0.312  0.17 

0.15 

0.16 

0.15 

0.18 

0.58 

0.48 

0.52 

0.53 

0.62 

Total  0.99 

0.98 

0.92 

0.98 

0.92 

0.66 

0.59 

0.59 

0.61 

0.61 

4.2.  Effect  of  the  kinetic  mechanism 

As  mentioned  previously,  since  in  flameless  combustion 
regime  the  turbulence— chemistry  interaction  as  well  as  the  ki¬ 
netic  mechanism  play  a  fundamental  role,  the  EDC  model  coupled 
with  four  different  kinetic  schemes  (JL,  KEE58,  GRI  2.11  and  GRI 
3.0)  was  considered.  For  all  cases  the  turbulence  model  was  the 
MKE. 

The  observation  of  axial  and  radial  velocity  profiles  indicated 
that  the  kinetic  scheme  does  not  have  a  great  influence  the  flow 
field  predictions. 


The  predictions  of  the  temperature  field,  however,  are  strongly 
affected  by  the  choice  of  the  kinetic  scheme.  This  is  well  evident  from 
Fig.  6(b),  in  which  the  measured  temperatures  along  the  axis  of  the 
chamber  are  compared  with  those  predicted  with  the  different  ki¬ 
netic  schemes.  The  GRI  2.11  is  the  mechanism  that  performs  better 
(average  deviation  1.8%),  although  also  the  KEE58  gives  satisfactory 
results  (average  deviation  2%),  confirming  the  necessity  of  taking 
into  account  at  least  skeletal  kinetic  mechanisms.  The  JL  scheme 
over-predicts  the  temperature  along  the  chamber  axis  (average  de¬ 
viation  7%),  whereas  the  GRI  3.0  under-predicts  it  (average  deviation 
8%).  Even  though  the  GRI  3.0  is  a  detailed  kinetic  mechanism,  it  fails 
to  predict  flameless  conditions,  as  also  reported  by  Sabia  et  al.  [56], 

4.3.  Effect  of  wall  boundary  condition 

The  heat  exchange  at  the  wall  has  a  major  effect  on  the  tem¬ 
perature  profile  inside  the  furnace.  During  the  experimental 
campaign,  the  wall  temperature  at  a  distance  z  =  0.112  m  was 
measured  as  1313  K  (1040  °C)  [29].  Being  flameless  combustion 
characterized  by  a  uniform  temperature  field,  it  was  considered 
acceptable  to  carry  out  numerical  simulations  with  walls  at  1313  K. 
In  order  to  verify  the  goodness  of  such  hypothesis,  also  a  convective 


Radial  distance  [m]  Radial  distance  [m] 


Radial  distance  [m]  Radial  distance  [m] 

Fig.  4.  Comparison  between  experimental  radial  profiles  of  the  axial  component  of  velocity  fluctuations  and  those  predicted  with  different  turbulence  models,  at  different  axial 
coordinates.  EDC  combustion  model,  KEE58  kinetic  mechanism. 
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(a)  (b) 

Fig.  5.  Parity  plot  of  measured  and  predicted  values  of  axial  (a)  mean  velocity  and  (b)  velocity  fluctuation,  for  different  turbulence  models.  EDC  combustion  model,  KEE58  kinetic 
mechanism. 


Table  4 

Coefficient  of  determination  (R2)  for  axial  velocity  and  axial  component  of  velocity 
fluctuations  with  different  turbulence  models.  EDC,  KEE58. 


MKE 

SKE 

RNG 

RKE 

RSM 

Axial  velocity 

Velocity  fluctuations  (RMS) 

0.826 

0.787 

0.893 

0.869 

0.880 

0.645 

0.846 

0.752 

0.876 

0.775 

heat  transfer  at  the  wall  was  considered.  First  a  mean  heat  transfer 
coefficient  was  calculated,  on  the  basis  of  a  global  heat  balance  on 
the  furnace.  The  mean  heat  transfer  coefficient  is  7  W/m2  K, 
compatible  with  a  natural  convection  process. 

Fig.  7  shows  the  temperature  profile  at  different  locations  along 
the  furnace  wall  using  the  heat  transfer  coefficient  boundary  con¬ 
dition.  It  can  be  observed  that  the  temperature  ranges  from  1320  K 
to  1380  K,  indicating  that  a  temperature  gradient  exist  and  the 
temperature  is  not  strictly  constant  at  the  wall.  However,  the 
relative  temperature  variation  between  the  bottom  and  the  top  of 
the  furnace  is  below  5%,  thus  indicating  that  the  isothermal  hy¬ 
pothesis  for  the  wall  is  reasonable. 

This  is  further  confirmed  by  the  analysis  of  Fig.  8,  which  compares 
the  predictions  of  the  temperature  along  the  axis  of  the  furnace 
when  the  two  different  boundary  conditions  are  considered.  It  is  not 
possible  to  distinguish  between  the  results  obtained  in  the  two  cases, 
indicating  that  the  type  of  boundary  condition  used  at  the  wall  does 
not  significantly  affect  the  predictions  inside  the  furnace. 


X  [m] 

(a) 


4.4.  NO  emissions 

As  shown  in  the  previous  section,  using  different  combustion/ 
kinetic  models  leads  to  different  temperature  (and  chemical  spe¬ 
cies)  distributions  in  the  furnace.  As  a  consequence  NO  emissions 
are  predicted  in  a  different  way. 

The  experimental  value  for  the  NO  emissions  of  the  furnace 
(provided  by  Plessing  et  al.  [28])  is  10  ppmv  on  dry  basis.  NO  (on 
dry  basis)  predicted  using  different  models  are  reported  in 
Table  5.  The  CR1  2.11  and  the  KEE58  schemes  perform  similarly 
and  yield  to  the  best  agreement  with  experimental  NO  data, 
predicted  NO  being  respectively  8.1  and  8.7  ppmv  and  experi¬ 
mental  ones  of  10  ppmv.  Conversely,  the  JL  scheme  provides  NO 
emissions  of  39.9  ppmv,  thus  far  from  the  experimental  evi¬ 
dences,  whereas  the  GR1  3.0  under-predicts  the  NO  emissions  by 
about  50%  (5.3  ppmv). 

Given  the  low  temperatures,  thermal  NO  formation  does  not 
play  a  major  role.  The  contribution  of  the  different  pathways  to  the 
overall  NO  emissions  is  shown  in  Fig.  9.  It  can  be  noticed  the 
dominant  role  of  the  N20  mechanism.  In  particular  for  the  cases  of 
GRI  2.11,  KEE58  and  GRI  3.0  models  the  N20  mechanism  is 
responsible  for  more  than  80%  of  the  NO  production.  With  the  JL 
scheme,  however  temperatures  are  higher  so  that  the  thermal 
mechanism  dominates  NO  formation  (99%  of  the  NO  production). 
This  is  due  to  the  completely  different  prediction  on  the  tempera¬ 
ture  field  inside  the  furnace. 


X  [m] 

(b) 


Fig.  6.  Comparison  between  experimental  temperature  profile  along  the  combustion  chamber  axis  and  those  predicted  with  (a)  different  turbulence  models  and  (b)  different 
kinetic  mechanisms.  MKE  turbulence  model,  EDC  combustion  model. 
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X  [m] 

Fig.  7.  Temperature  profiles  at  different  locations  along  the  furnace  wall,  predicted 
with  a  convective  heat  transfer  wall  boundary  condition.  MKE  turbulence  model,  EDC 
combustion  model,  KEE58  kinetic  mechanism. 

5.  Uncertainty  analysis 

Validation  cannot  be  carried  out  without  explicitly  accounting 
for  the  uncertainties  present  in  both  the  measurements  and  the 
computation.  As  far  as  computation  is  regarded,  the  uncertainties 
are  associated  to  the  choice  of  the  physical  model  and  to  the 
specification  of  the  input  parameters  required  for  performing  the 
analysis.  The  objective  of  the  present  section  is  to  show  how  sen¬ 
sitive  the  numerical  simulations  are  with  respect  to  scenario  pa¬ 
rameters  (boundary  conditions)  and  model  parameters.  This 
appears  crucial  to  identify  the  main  sources  of  uncertainties  and 
orient,  therefore,  modelling  (and  experimental)  improvements. 

The  specification  of  the  boundary  conditions  is  a  key  issue  in 
numerical  simulation,  and  typically  only  limited  information  is 
available  from  corresponding  experiments  and  observations.  In  the 
case  of  the  furnace  under  investigation,  the  major  uncertainty  is 
related  to  the  wall  boundary  condition.  During  the  experimental 
campaign,  the  wall  temperature  at  a  distance  z  =  0.112  m  was 
measured  as  1313  K  (1040  °C)  [29],  Being  flameless  combustion 
characterized  by  a  uniform  temperature  field,  it  was  considered 
acceptable  to  carry  out  numerical  simulations  with  walls  at  1313  K, 
as  discussed  in  Section  3.1. 


X  [m] 

Fig.  8.  Comparison  between  experimental  temperature  profile  along  the  combustion 
chamber  axis  and  those  predicted  with  different  wall  boundary  conditions.  MKE  tur¬ 
bulence  model,  EDC  combustion  model,  KEE58  kinetic  mechanism. 


Table  5 

Comparison  between  experimental  NO  emissions  [ppm„]  and  those  predicted  using 
different  kinetic  schemes. 


Exp. 

JL 

KEE58 

GRI2.11 

GRI3.0 

10 

39.9 

8.7 

8.1 

5.3 

We  can  now  assess  the  sensitivity  of  the  numerical  results  to 
such  a  choice  by  modifying  the  wall  temperature  within  a 
reasonable  range  of  variability.  Following  the  recent  analysis  pro¬ 
vided  by  Parente  et  al.  [55  ,  who  analyzed  the  possible  uncertainty 
sources  for  in-furnace  temperature  measurements,  an  uncertainty 
of  ±40  K  is  considered  here.  Therefore,  two  temperature  values 
were  used  to  specify  the  furnace  wall  temperatures  beside  the  base 
case  discussed  above:  1273  I<  and  1353  K.  All  the  simulations  were 
run  using  the  modified  k—e  turbulence  model  and  the  EDC  com¬ 
bustion  model  coupled  with  the  KEE58  chemical  mechanism,  since 
it  allows  having  good  results  and  saving  computational  time  with 
respect  to  GRI  2.11. 

Fig.  10(a)  shows  the  axial  temperature  profiles  for  these  three 
cases.  It  is  possible  to  notice  that  the  wall  temperature  has  a  major 
effect  on  the  temperature  inside  the  chamber.  In  this  case  the  best 
agreement  with  the  experimental  data  is  obtained  when  the  value 
measured  in  the  experimental  campaign  was  considered. 

In  particular,  when  a  lower  temperature  is  considered,  the 
maximum  temperature  inside  the  chamber  decreases  of  125  I<  with 
respect  to  the  baseline  case  (Tw  =  1313  K),  leading  to  an  average 
deviation  of  9%.  On  the  other  hand  increasing  the  wall  temperature 
leads  to  an  average  deviation  of  7%  and  an  increase  of  the  maximum 
temperature  of  80  K.  Therefore,  the  propagation  of  the  wall  tem¬ 
perature  uncertainty  in  the  simulation  results  is  slightly  more  than 
linear. 

Following  our  discussion  in  Paragraph  4,  the  uncertainty  related 
to  the  heat  transfer  coefficient,  h,  at  the  wall  was  also  considered.  A 
variation  of  10%  for  h  was  considered  therefore  two  additional 
simulations  were  run  with  h  =  6.3  W/m2  K  and  h  =  7.7  W/m2  K, 
respectively.  Results  are  provided  in  Fig.  10(b).  When  a  lower  co¬ 
efficient  is  used,  the  temperature  inside  the  furnace  increases  of 
90  K,  while  with  the  higher  coefficient  the  temperature  decreases 
of  60  K,  thus  showing  that  the  propagation  of  the  uncertainty  is 
slightly  less  than  linear.  This  is  an  interesting  result,  showing  that  a 
convective  heat  transfer  boundary  conditions  might  be  preferable 
over  a  fixed  temperature  one,  since  an  error  in  the  specification 
boundary  condition  would  have  a  less  significant  effect  on  the 


Fig.  9.  Contribution  of  different  formation  routes  to  the  total  NO  emissions  (dry  basis), 
predicted  with  different  kinetic  mechanisms.  MKE  turbulence  model,  EDC  combustion 
model. 
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Fig.  10.  Uncertainty  on  axial  temperature  predictions,  associated  with  (a)  the  wall  temperature  boundary  conditions,  (b)  the  heat  transfer  coefficient,  (c)  the  turbulence  model  and 
(d)  the  kinetic  mechanism.  EDC  combustion  model,  KEE58  kinetic  mechanism. 


results.  The  mean  wall  temperature  is  1430  K  with  h  =  6.3  W/m2  K 
and  1304  I<  with  h  =  7.7  W/m2  K. 

It  is  also  possible  to  compare  the  effect  of  the  uncertainty  in  the 
boundary  conditions  to  the  one  related  to  the  choice  of  the  physical 
models  (i.e.  turbulence,  combustion  and  kinetic  model). 

Fig.  10(c)  and  (d)  shows  the  variability  in  the  results  related  to 
the  variation  of  the  physical  model  and  the  kinetic  mechanism, 
respectively.  The  uncertainties  related  to  the  boundary  conditions 
and  to  the  physical  model  are  of  the  same  order  of  magnitude.  It  is 
possible  to  notice  a  difference  of  120  I<  on  the  maximum  temper¬ 
ature  inside  the  chamber  as  a  result  of  the  change  in  the  physical 
model.  The  average  deviation  is  2%  for  MKE  and  8%  for  SKE. 

A  larger  uncertainty  is  associated  with  the  choice  of  the  kinetic 
mechanism.  Switching  from  JL  to  GRI  3.0  leads  to  a  difference  on 
the  maximum  temperature  of  about  180  K.  The  relative  average 
deviation  is  1.8%  for  GRI  2.11,  2%  for  KEE58,  7%  for  JL  and  8%  for  GRI 
3.0.  This  confirms  the  importance  of  the  choice  of  the  kinetic 
mechanism  and  in  the  specific  case  of  flameless  furnaces,  the  ne¬ 
cessity  of  taking  into  account  detailed  chemistry. 


First  a  theoretic  estimation  of  I<r  is  performed,  following  Cardoso 
[57],  who  relates  recirculation  rate,  kRt  to  the  chamber-to-nozzle(s) 
areas  by  means  of  the  parameter  a: 

kR  =  0.30aVl  59  (4) 


(5) 


where  Ac  is  the  section  of  the  combustion  chamber  and  An  the 
section  of  the  inlet  nozzles.  Using  Equation  (5),  a  value  of  9.7  is 
obtained  for  kR,  in  agreement  with  the  value  provided  by  provided 
by  Plessing  et  al.  [28] 

The  value  of  I<r  estimated  using  Equations  (4)  and  (5)  can  be 
verified  by  post-processing  the  CFD  results.  Starting  from  the  flow 
streamlines  inside  the  furnace,  it  is  possible  to  locate  the  plane  on 
which  the  main  recirculation  takes  place.  The  value  of  kR  is  then 
calculated  from  the  following  equation,  reported  in  Ref.  [1[: 


6.  Operating  characteristics  of  the  furnace 


1<R  = 


me 
rrif  +  ma 


(6) 


Once  the  computational  model  is  chosen  and  the  uncertainties 
related  to  both  models  and  boundary  conditions  are  assessed,  it  is 
possible  to  evaluate  some  of  the  operating  characteristics  of  the 
furnace,  such  as  the  recirculation  rate  kR  and  the  location  of  the 
reaction  zone,  in  order  to  show  that  some  of  the  theoretical  char¬ 
acteristics  of  the  flameless  regime  are  actually  present  in  the  case 
under  investigation. 

In  flameless  combustion  regime,  the  recirculation  rate  is  a  key 
parameter  to  quantify  the  amount  of  exhaust  gases  recirculated  [1  ]. 
The  recirculation  rate  is  defined  as  the  ratio  between  the  net  mass 
flow  rate  of  recirculated  flue  gas  and  the  sum  of  the  fuel  and  air 
mass  flow  rates. 


where  me  is  the  net  mass  flow  rate  of  recirculated  flue  gas,  whereas 
rrif  and  ma  are  the  fuel  and  air  mass  flow  rates,  respectively.  All 
those  mass  flow  rates  are  calculating  following  Equation  (7): 

m  =  J  pvdS  (7) 

s 

The  value  provided  by  Equations  (5)  and  (6)  is  kR  =  10.7,  in 
excellent  agreement  with  the  theoretical  value  calculated  above. 

In  the  conditions  of  the  traditional  combustion  regime  it  is  not 
possible  to  achieve  flammable  mixtures  of  hydrocarbon  and  air  for 
values  of  kR  >  0.5  without  extinction  occurring,  due  to  the  lower 
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oxygen  concentration  and  higher  inert  species  in  the  mixture. 
However,  the  temperature  and  mixing  ensured  in  the  furnace  allow 
the  fuel  to  react  in  a  steady  and  stable  form  even  for  high  values  of 
the  recirculation  rate;  this  is  the  principle  behind  flameless  com¬ 
bustion  [1,2], 

As  far  as  the  reaction  zone  is  concerned,  it  is  possible  to  locate  it 
considering  the  distribution  of  two  chemical  species:  the  radical 
OH  and  the  formaldehyde  (CH2O).  Fig.  11  shows  the  distribution  of 
the  mass  fraction  of  OH  and  CH2O  across  the  furnace  axis.  It  is  clear 
that  the  reaction  zone  is  lifted,  meaning  that  the  reaction  is  taking 
place  in  diluted  condition  and  not  as  soon  as  the  fuel  and  air  jets 
mix. 

In  order  to  verify  if  the  reaction  is  actually  taking  place  in  diluted 
condition  the  distribution  of  the  mass  fraction  of  OH  versus  the 
mass  fraction  of  O2  is  analyzed.  Fig.  12  compares  such  distribution 
for  three  different  flames: 

-  Flame  C  of  the  Sandia  Laboratory  58],  which  is  a  purely  diffusive 

flame; 

-  JHC  burner  [21  ,  which  emulates  flameless  conditions; 

-  Furnace  under  investigation. 

As  far  as  the  first  two  systems  are  concerned,  the  experimental 
data  were  used,  whereas  for  the  furnace  under  investigation  the 
data  obtained  by  the  numerical  model  validated  in  the  previous 
sections  are  used. 

It  is  possible  to  observe  that  for  flame  C  (red  stars,  in  the  web 
version)  the  reaction  takes  place  for  all  the  possible  concentration 
of  O2,  whereas  for  the  JHC  (blue  crosses,  in  the  web  version)  the 
reaction  takes  place  only  in  very  diluted  conditions.  The  furnace 
under  investigation  is  a  flameless  system.  However,  a  flameless-like 
behaviour  can  be  only  observed  with  an  appropriate  choice  of 
combustion  model  and  kinetic  mechanism.  To  illustrate  this,  the 
mass  fraction  obtained  with  the  KEE  mechanism  is  plotted  versus 
the  mass  fraction  of  O2  (black  circles).  It  can  be  observed  that  the 
combination  EDC  with  KEE  is  able  to  predict  that  reactions  mostly 
take  place  in  diluted  condition  (mass  fraction  of  02  below  0.06).  A 
relatively  small  region  where  the  reactions  take  place  in  presence  of 
higher  concentration  of  O2  is  also  identified:  the  latter  is  probably 
associated  to  the  ignition  process,  which  is  likely  controlled  by 
auto-ignition  and  takes  place  in  diffusive  conditions,  as  supported 
by  recent  investigations  on  the  flameless  combustion  regime  [59], 


X  [m] 


Fig.  11.  OH  and  CH20  mass  fraction  distribution  along  the  furnace.  MKE  turbulence 
model,  EDC  combustion  model,  KEE58  kinetic  mechanism. 


Fig.  12.  Distribution  of  the  mass  fraction  of  OH  versus  the  mass  fraction  of  02  for  Flame 
C,  JHC  and  the  furnace  under  investigation. 


Therefore,  we  can  conclude  that  the  approach  described  in  the 
previous  sections  lead  to  a  comprehensive  numerical  model, 
which  well  describes  the  actual  combustion  system  as  indicated 
by  the  quantitative  agreement  between  simulations  and  experi¬ 
mental  observations.  It  should  be  stressed  that  the  model 
development  was  not  based  on  qualitative  assessment  and  trial 
and  error  procedures,  but  it  has  been  guided  by  a  rigorous 
analysis  of  the  possible  sources  of  modelling  and  scenario 
uncertainties. 

7.  Concluding  remarks 

A  numerical  investigation  through  computational  fluid  dy¬ 
namics  of  a  semi-industrial  furnace  operating  in  flameless  com¬ 
bustion  mode  has  been  presented.  In  particular  the  role  of  closure 
sub-models  in  flow  field,  temperature  and  NO  emissions  was 
studied. 

Among  the  different  turbulence  models  that  were  tested,  the 
modified  k—e  model  provides  the  best  agreement  with  the  exper¬ 
imental  data,  especially  as  far  as  temperature  is  concerned.  Rey¬ 
nolds  stress  model  leads  to  a  smaller  deviation  from  the  measured 
value  when  radial  velocity  and  fluctuation  velocity  are  concerned. 

The  Eddy  Dissipation  Concept  (EDC),  coupled  with  four  different 
kinetic  schemes  (JL,  KEE58,  GRI  2.11  and  GRI  3.0)  was  considered. 
The  GRI  2.11  is  the  mechanism  that  performs  better,  although  also 
the  KEE58  gives  satisfactory  results,  confirming  the  necessity  of 
taking  into  account  at  least  skeletal  kinetic  mechanisms. 

As  far  as  NO  emissions  are  concerned,  similar  results  were 
found.  The  GRI  2.11  and  the  KEE58  perform  similarly  and  yield  to 
the  best  agreement  with  experimental  NO  data. 

The  contribution  of  different  formation  routes  was  also  exam¬ 
ined.  Given  the  low  temperatures  that  characterize  such  combus¬ 
tion  regime,  the  N2O  intermediate  mechanism  for  formation  of  NO 
plays  a  major  role,  while  thermal  NO  mechanism  is  not  as  relevant 
as  in  traditional  combustion  regime. 

Finally  an  assessment  of  the  uncertainty  related  to  the  choice  of 
boundary  conditions  was  performed.  First,  the  effect  of  two 
different  thermal  boundary  conditions  was  evaluated  and  then  the 
uncertainty  related  to  the  value  of  the  wall  temperature  and  the 
value  of  the  heat  transfer  coefficient  was  assessed.  It  was  found  that 
the  wall  temperature  has  a  major  effect  on  the  temperature  inside 
the  chamber  and  this  uncertainty  is  of  the  same  order  of  magnitude 
as  the  uncertainty  associated  with  the  choice  of  the  physical 
models.  A  larger  uncertainty  is  associated  to  the  kinetic  mecha¬ 
nism,  confirming  that  it  plays  a  major  role  in  the  modelling  of 
flameless  furnaces. 
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Finally  the  operation  characteristics  of  the  furnace  are  evalu¬ 
ated.  The  recirculation  rate  I<r  is  estimated  at  10.7  and  it  was  found 
from  the  analysis  of  the  radical  OH  distribution  that  the  reaction 
takes  place  mostly  in  very  diluted  conditions. 

The  model  obtained  in  the  present  work  well  agrees  with  the 
experimental  data  and  it  is  the  result  of  a  rigorous  analysis  of  the 
possible  sources  of  modelling  and  scenario  uncertainties. 
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